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Abstract 
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QCD with a finite baryon chemical potential, despite its importance, is not well understood 
CN because the standard lattice QCD simulation is not applicable due to the sign problem. Although 

QCD-like theories which do not suffer from the sign problem have been studied intensively, relation 
to QCD with a finite baryon chemical potential was not clear. This paper introduces \aige-Nc 
'"P' equivalences between QCD and various QCD-like theories. These equivalences lead us to a unified 

^ viewpoint for QCD with baryon and isospin chemical potentials, SO{2Nc) and Sp{2Nc) gauge 

(— I theories, QCD with adjoint matters and two-color QCD. In particular QCD with the baryon 

chemical potential is equivalent to its phase quenched version in a certain parameter region, 
1— I which is relevant for heavy ion collision experiments. This result is useful especially for finding 

<^ or excluding a possible QCD critical point, through the phase quenched simulation of the SU{S) 

^ QCD. 
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1 Introduction 



Consider QCD at a finite baryon chemical potential, 

1 _ 

^ /=1 

where the gauge group is SU{3), Nf is the number of flavors, ipf are quarks of mass in 
the fundamental representation, and /x is the quark chemical potential which is related to the 
baryon chemical potential fiB as fiB = 3/x. Properties of this theory have long been a subject of 
intense interest. (For a review, see [Ij.) Apart from its intrinsic theoretical appeal, this subject 
is important in astrophysics, especially in the study of neutron stars. The behavior of QCD at 
asymptotically large //_b is well understood theoretically due to the asymptotic freedom, and QCD 
becomes a color superconductor as //s — )■ oo [2]. 

At more phenomenologically realistic densities, QCD is strongly coupled, and thus not amenable 
to controlled analytic treatment. Although lattice Monte Carlo is very useful at = 0, how- 
ever, it runs into trouble at fis 7^ due to the fermion sign problem - the fermion determinant 
ri/iii "^^^ {l^^fi + rnj + ^7^) becomes complex, rendering importance sampling exponentially dif- 
ficult. 

In order to circumvent this difficulty people have studied gauge theories which do not suffer 
from the sign problem at finite density. Consider QCD and QCD-like theorieaj of the form 




where G is the gauge group e.g. SU{Nc), ^/ is a generic quark chemical potential, and fermions 
■0*-' are not necessarily in the fundamental representation. The main examples are QCD with an 
isospin chemical potential /// (i.e. Nf = 2, /ii = —^2 = fJ^i/^) and degenerate mass mi = m2 [3], 
two-color QCD of even number of flavors and degenerate mass QCD with adjoint fermions [1], 
and SO{2Nc) and Sp{2Nc) Yang-Mills theories [5l|6]. However, while interesting, these theories 
have many qualitative differences from Nc = 3 QCD, such as e.g. explicitly broken flavor symmetry 
in the first case. Therefore it is important to understand what we can learn from these theories, 
or in other words, in what sense they are similar to real QCD with the baryon chemical potential. 
In [3 [B], an answer to this question has been given. The statements are 

• SO{2Nc) YM with hb, Sp{2Nc) YM with ^lB and SU{Nc) QCD with in are large-iVc 
equivalent both in the 't Hooft limit {Nf fixed) and the Veneziano limit (Nf/Nc fixed), 
everywhere in the T-^ plane. (Fermions are in the fundamental (vector) representation.) 

• SO{2Nc) YM with fiB, Sp{2Nc) YM with ^b, SU{Nc) QCD with fij and SU{Nc) QCD with 
^B are large- A/'c equivalent in the 't Hooft limit, outside the BEC/BCS crossover region of 
the former three theories. (Fermions are in the fundamental representation.) 

^ In this paper we call SU{Nc) Yang-Mills with Nf fundamental fermions 'QCD'. SU{Nc) Yang-Mills with 
fermions in other representations and SO{2Nc) / Sp{2Nc) theories are referred to 'QCD-like theories'. 
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• More generally, SO{2Nc), Sp{2Nc) and SU (Nc) theories with fermion mass mi, • • • ,mNf and 
chemical potential /ii, • • • jfJ-Nf are equivalent. The signs of the quark chemical potentials 
can be flipped without spoiling the equivalence. (Fig. [2]) 

• SO{2Nc) YM with the Nf complex adjoint fermions and fiB, SU{Nc) YM with the Nj 
complex adjoint fermions and fis, and SU{Nc) YM with the 2Nf complex anti-symmetric 
fermions and /U/ are large- A'^c equivalent everywhere in the T-fj, plane. 

These statements have been derived by using string-inspired large- A'^c techniques [71 [HI l9l [TO]. 
As shown in O |6], there are orbifold projections relating S0{2Nc) and Sp{2Nc) theories with 
Hb, QCD with hb and QCD with iij (Fig. [T]). At large-A^ci the orbifold equivalence guarantees 
these theories are equivalent in the sense a class of correlation functions (e.g. the expectation 
value of the chiral condensate and vr^ correlation functions) and the phase diagrams determined 
by such quantities coincide, as long as the projection symmetry is not broken spontaneously |10j . 
A similar argument shows QCD with adjoint fermions and /i_B is equivalent to QCD with fermions 
of two-index antisymmetric representation, which is a so-called the Corrigan-Ramond large-iVc 
limit, with m (Fig. [s]). In order for these equivalence to hold, orbifolding symmetries must 
not be broken spontaneously. This requirement is always satisfied for the equivalences between 
S0{2Nc) YM with /i^, Sp{2Nc) YM with fiB and QCD with isospin chemical potential. For the 
equivalences between these three theories and QCD with the baryonic chemical potential, 'outside 
the BEC/BCS crossover region' is required for the symmetry realization. This region ('outside 
the BEC/BCS crossover region') is relevant for the search for the 'QCD critical point', which 
attract intense interest over the decade. Our answer to the problem is strikingly simple - by using 
sign-free theories one can answer to the question. In the case of the two-flavor theory, QCD with 
m is nothing but the phase-quenched version of QCD with fiB- Therefore, the sign problem is 
merely an illusion, up to the 1/Nc correction. Furthermore, for gluonic observables, the leading 
1/Nc corrections to the large- A^c limit which carries the information of the chemical potential are 
the same in these theories. We also argue that it might be possible to study the large- //^ region, 
where the equivalence fails because of the symmetry breakdown, by adding certain deformations 
to the SO{2Nc) theory while leaving the QCD with hb untouched. 

This paper is organized as follows. In § [2] we show the absence of the sign in S0{2Nc) theory 
(and also in Sp{2Nc) theory when the mass and chemical potential are degenerated). In § |3| we 
review the orbifold equivalence and provide a proof to all order in perturbatioij^ As discussed in 
|10) . projection symmetries must be unbroken in order for the equivalences to hold. In order to 
see the fate of the symmetries, we discuss the phase diagrams of S0{2Nc), Sp{2Nc) and SU{Nc) 
theories in § [3.2.3 In § [4] we generalize the equivalences to analytically solvable toy models and 



confirm the equivalence explicitly. Results in this section strongly suggests the validity of the 
equivalences at nonperturbative level. § [5] is devoted for the conclusion and outlook. 

This paper has been prepared for Seitaro Nakamura prize competition, and is based on papers 
[5] and [6]. We also provide a few new results. 



^ We do not show the proof the equivalence between theories with adjoint and antisymmetric fermions. The 
proof for the theories with fundamental fermion can easily be modified for this case. The projection is the same as 
the one used in The equivalence between two SU{Nc) theories (when the chemical potential is zero) has been 
found in [H]. 
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S0(2Nc) YM^ ^Sp(2Nc) YM 

with |iB with |iB 




SU(Nc)QCD^ ^SU(Nc)QCD 
with |iB ' * with |i, 



Figure 1: A web of equivalences. Arrows with solid lines represents equivalences through orbifold 
projections. Arrows with dashed lines are the 'parent-parent' and 'daughter-daughter' equivalences 
which arise as combinations of the orbifold equivalences. 



S0(2Nc) YM^ ^Sp(2Nc) YM 

with [i2,... with |ii, |i2,... 




SU(Nc)QCD 
with |i2,... 



Figure 2: More general version of the equivalences. Values of the quark chemical potentials can 
be different. 
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S0(2Nc), Nf adjoint 
with |iB 




SU(Nc), SU(Nc), 

Nf adjoint 4 ► 2NfAS 

with [Xb with |i, 

Figure 3: Equivalences in the Corrigan-Ramond limit. SU (A^c) YM with anti-symmetric fermions 
can be regarded as a special kind of large- A/'c limit of three-color QCD (the Corrigan-Ramond 
limit), because anti-symmetric and fundamental representations are equivalent at Nc = 3. Un- 
fortunately, SU{Nc) YM with anti-symmetric fermions and /ib cannot be incorporated in these 
equivalences. 

2 Absence of the sign problem in SO{2Nc) and Sp{2Nc) theories 

In this section we prove the absence of the sign problem in SO{2Nc) and Sp{2Nc) theories. 
Let us start with SU{Nc) theory. In four-dimensional Euclidean space, the gamma matrices 
can be taken to be Hermitian, = 7/*- The covariant derivative = 5^ + , where 

is an Nc x Nc Hermitian matrix, is skew- hermitian, (D^^)^ = —D^. Therefore, 7//-D/i is also 
skew- Hermitian, 

(7'^0^ = -7'^<- (3) 

Hence the eigenvalues of 7^D^^ are pure imaginary. Furthermore the eigenvalues appear in pairs 
±iX, where A is real, because of the chiral symmetry 75 (7'*!?^^) 75 = —7'*!?^^. When the mass 
irif is turned on, the eigenvalues are shifted to itiA + mj. Hence the eigenvalues appear with their 
complex conjugate (as long as m/ is real) and hence the determinant, which is the product of the 
eigenvalues, is real positive. 

Once the chemical potential 11 is turned on, this property is lost because 

irDf" + /X74) ^ = -I'^Df' + IJilA, (4) 
and hence the determinant is complex in general. Note however that 

[det (j^'D^ + m + M7^) ] * = det {j^'D^ + m - /i^^) . (5) 
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For this reason QCD with the isospin chemical potential and degenerate mass is sign-free; the 
determinant is 



det {-i^D^ + m + fi-/^) x det {"t^D^, + m - fi'y^) = |det {-f'^D^ + m + ^-f^) f >0. (6) 

In SO{2Nc) theory, the sign is absent thanks to an additional anti-unitary symmetry. The 
crucial point is -D^^ is real in the coordinate basis; both 9^ and ij4^^ are real antisymmetric. 
Therefore 

(C75) {l^'D^ + m + (C75)-' = (7^^M + ^ + /^t')* , (7) 

where C is the charge conjugation matrix satisfying C^^C~^ = — = ~7^- It guarantees the 
pair structure of the eigenvalues (A, A*); if v is an eigenvector with an eigenvalue A, {C^^)~^v* is 
another eigenvector with an eigenvalue A*, and furthermore, they are linearly independent even 
when A is real [S]. 

In Sp{2Nc) theory, a similar relation holds, 

( J,C775) {^D^ + m + /X74) ( JeC75)-' = {l^D^ + m + fi^^ , (8) 



where Jc = i(T2 ® Iat^ (see § 3.1.2). Unlikely to the case of SO{2Nc), this relation does not give 
the pair structure when A is real, and hence only the reality, not the positivity, is guaranteed. 
However as long as the mass and chemical potential are degenerate (mi = m2, ms = uia, • • • , 
IJ-i = ^J'2, ^J'3 = fJ-i, • • •) the determinant is real and positive. 



3 Orbifold equivalence 

3.1 Pure Yang-Mills theory 

3.1.1 SO{2Nc) to SU{Nc) 

The notion of the orbifold equivalence came from the string theory [7, 8j. Soon it has been 
proven by using only field theory techniques [9l |10] , without referring to the string theory. As a 
simple example, let us consider the equivalence between SO{2Nc) and SU{Nc) pure Yang-Mills 
theories. To perform an orbifold projection, one identifies a discrete subgroup of the symmetry 
group of the 'parent' theory, which is the S0{2Nc) theory in this case, and require the fields to 
be invariant under the discrete symmetry. This gives a 'daughter' theory, which is SU{Nc) YM. 
The orbifold projection uses a Z2 subgroup of the S0{2Nc) gauge symmetry. 

Let us take Jc € S0{2Nc) to be Jc = ia2 Iat^, which generates a Z4 subgroup of S0{2Nc). 
Here Iat is an x identity matrix. We require the gauge field to be invariant under 

-Afi ^ JcAfxJf, , (9) 

which generates a Z2 subgroup of S0{2Nc). A^^ can be written in Nc x Ac blocks as 

'~\C^ + D'', A^-B^)^ ^^0^ 
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where fields with an 'A' ('5") superscript are anti-symmetric (symmetric) matrices. Under the 
Z2 symmetry, A'^,D^ are even while B^,C^ are odd, and hence the orbifold projection sets 
= C^ = 0; the 'daughter' field is 

^r=M j^'s ^/ ■ (11) 



By using a unitary matrix 



it can be rewritten as 



p _ JL l^e ^liv, I 

V2 V 1^. -^1^. ' ' ^ ^ 



where = + is a U{Nc) gauge field. However, the difference between U{Nc) and SU{Nc) 
is a correction and can be neglected at large- A^'c. The gauge part of the action after the 

orbifold projection is thus simply 

^gaugcproj ^ ^_ J- ^^M^^ (14) 

where is the SU{Nc) field strength. Let us identify it with the Lagrangian of the daughter 
theory times two, 

£so ^ 2£su, (15) 

or equivalently let us take g^jj = g^Q, where gsu is the gauge coupling constant of the SU{Nc) 
theory. This factor two is necessary in order for the ground state energies, which are proportional 
to the degrees of freedom, to match. Then expectation values of the gauge-invariant operators in 
parent agree with the expectation values of the daughter theory, which is obtained by 

replacing A^, with O'^^^iAf,] = 

As a pedagogical demonstration, consider a planar diagram in Fig. |4j In order to obtain the 
SU{Nc) diagram, we insert the projector V defined by 



Ar = nAf) = \ jz J^Al^J-- = \ [Af + JcAl^J-') . (16) 

n=0 

to each propagator in the SO{2Nc) diagram. Then the only difference, if exists, comes from the 
contractions of color indices. This additional kinematic factor multiplied to the SU{Nc) diagram 
is 

i^-j ■ tr(J""i J"5) . tr(J-"2 j~"V"«) • tr(J~"3 J-"5 J-"6) . tr(J"V"2 J"^), (17) 



ni=0,l 



^ When one studies U{Nc) theory, it is difficult to control the (7(1) part and the effect of the chemical potential 
can disappear [1^ by a lattice artifact. In order to avoid that one should simulate SU{Nc) theory on the lattice. 
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jn2 

Figure 4: A vacuum planar diagram in the double-line notation. 

where J = —ia2 is a 2 x 2 matrix and the factor (1/2)^^ comes from the projectors with Np = 6, 
where Np is the number of propagators. Because J satisfies a simple relation 

Trr = unless J" = ila, (18) 

it is nonvanishing only when 

j-n,jn^jns ^ j-n^ j-n^ jn^ ^ ^-^^^ J-^V^^s J""" = ±13, J«i J"2 jns = ^J^, 

(19) 

or equivalently 

—rii + n4 + ns = even, — n2 — n4 + ng = even, — ns — ns — ng = even, ni + n2 + ns = even. 

(20) 

Not all constraints are independent; actually the last one follows from the others, and hence, there 
are N^ — 1 = 3 independent constraints, where N^ = 4 is the number of index loops. By summing 
over m for all Np = 6 propagators with A'^^ — 1 = 3 constraints, one obtains a factor of 2^^^. 
Another factor 2^ comes from the traces over color indices. Therefore, the total factor is 

2-6 . 26-3 . 24 = 2. (21) 

Generally, for given planar vacuum diagrams with Np propagators and Ni loops, the projectors 
give a factor of (1/2)^^, the summation over rii under the Nl — 1 constraints gives 2^^^^^^^"'^^ 
and the trace gives 2^^. The total factor is always 2: 

2-iVp . 2^p-(7V^-i) . 2^L = 2. (22) 

This factor 2 reflects the fact that the number of degrees of freedom in the parent theory is 
twice larger than that in the daughter theory. Hence the vacuum energy per degree of freedom is 
equivalent between these theories. 

The counting does not apply for nonplanar diagrams. Indeed, one can easily check that the 



number of independent constraints is no longer Nl — 1, and the factor counted in (22) is generally 
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different from 2 [9j. This is why we need to take the large- A'c Umit to suppress the nonplanar 
diagrams. 

A few remarks are in order here. Firstly, in order for the orbifold equivalence to work, the 
projection symmetry must not be broken spontaneously. In this example, because the projection 
symmetry is embedded to the gauge transformation, it is not broken. In the following sections, 
we introduce matter fields and use flavor symmetries for the projection. Then the projection 
symmetry does break in a certain parameter region. Secondly, not all operators in two theories 
coincide. In the parent theory, only operators invariant under the projection symmetry is related 
to projected fields in the daughter. (In the above example this condition is not relevant because 
all the gauge invariant operators automatically satisfy this condition.) In the daughter, not all 
operators are obtained from the parent through the projection; in the present case, the projected 
operators are necessarily charge conjugation invariant. 



3.1.2 Sp{2Nc) to SU{Nc) 

For Sp{2Nc) gauge theory, the symplectic algebra Sp{2Nc) formed by 2Nc x 2Nc Hermitian 
matrices satisfying 

JcA^P + {A^Pfjc = 0, (23) 
can be written using Nc x Nc matrices as 

where the fields {B^, C^, and D^) are anti-symmetric (symmetric) matrices. After the 
projection, one obtains = = 0. This gives SU{Nc) gauge theory again. 



3.2 Introducing fundamental fermions 



3.2.1 Orbifold projection of fundamental fermions 

In this section, we introduce the orbifold projection for fundamental fermions l5j£j. 
Let us consider the effect of the orbifolding on ipf. By using oj = e™!'^ G [/(1)b, we define the 
projection by 

ipf = ujJcijjf. (25) 

Writing 

(I ) = (26) 
the action of the Z2 symmetry is just (V'j , V'/ ) ~^ (~^/ )■ '^^^ projection consists of setting 
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The action of the daughter theory is (after rescaUng the couphng constant) 



Nf 



(27) 



where J^f^^ is the field strength of the SU{Nc) gauge field = + iA^, ip^^ = and 
T^fi = + iAfj,. This is an SU{Nc) gauge theory with Nf flavors of fundamental Dirac fermions 
and the baryon chemical potential fiB = fJ'N'c- So the orbifold projection relates SO{2Nc) gauge 
theory to large Nc QCD. 

On the other hand, in order to obtain fermions at finite /x/ for even Nf, we use Jc G SO{2Nc) 
[or Je G Spi2Nc)] and Ji G 5C/(2)isospin[C SUiNf)] defined by 

Ji = -icr2'^lNf/2- (28) 
We choose the projection condition to be 

The flavor A^j-component fundamental fermion is decomposed into two (A^j/2)-component fields, 

V.so(Sp) = (^^ V',), (30) 



with i and j being the isospin indices. If we define = (V^j. =F ) / and ^± = (V'±i iV'4)/V2, 



the fermions ip± survive but ^± disappear after the projection (29). Because ip± couple to {A^)'^ 
and A^ respectively, the action of the daughter theory is expressed as 



Nf/2 



su 



tr(F,l^)^ + E E i^'D^ + - ± /^^') ^T±^ 



^9lv f=, ± 



(31) 



where V'S.^ = V^f- and = \/2ip'l. This theory has the isospin chemical potential /i/ = 2//. 



3.2.2 't Hooft limit vs Veneziano limit 



In the proof of the orbifold equivalence of the pure Yang-Mills theories shown in § 3.1, the 
conditions (20) are crucial. What happens when the fermions are introduced? Firstly note that 
two projections (25) and (29) are equivalent when the chemical potential is absent. Both are a Z4 



subgroup of the fiavor symmetry which mixes two Majorana flavors. 

Once the chemical potential is turned on, they are not equivalent. 

Therefore 



The flavor symmetry 
the proof can be repeated 
straightforwardly; the only difference is some color-index loops are replaced with flavor-index 



Ji used in (29) satisfles the condition similar to (18). 



loops. On the other hand, Z4 G [/ (1) b used in (25 ) does not satisfy such a condition. Note however 



that one of the conditions in (20) is not independent and follows from other ones. Therefore, as 



long as the number of flavor-index loop is one, the proof holds. Because the fermion loops are 



suppressed by the factor Nf/Nc, the equivalence through (25) holds in the 't Hooft large- A^c limit 



{Nf fixed) while the one through (29) holds also in the Veneziano limit (Nf/Nc fixed). 

The above argument has an implication for the 1/Nc correction. Let us consider QCD with 
fiB and that with ^j. In the 't Hooft large- A'^c limit, gluonic operators trivially agree because the 
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fermions are not dynamical. Let us consider finite- A'^ci say Nc = 3 and Nf = 2. Then the largest 
correction to the 't Hooft hmit comes from one-fermion-loop planar diagrams, which, as we have 
seen, do not distinguish fj,B and /x/. Therefore gluonic operators should behave similarly even 
quantitatively; the difference is at most {Nf /Nc)'^ (two-fermion-loop planar diagrams) or (l/N^) ■ 
(Nf/Nc) (one-fermion-loop nonplanar diagrams). In particular, the deconfinement temperatures, 
which is determined from the Polyakov loop, should be close. 



3.2.3 Symmetry realization and validity of the equivalence 

As we have seen so far, SU{Nc) QCD with fis/f^i, S0{2Ne) YM and Sp{2Nc) YM should be 
equivalent in the large- A'^c limit as long as the projection symmetries are not broken spontaneously. 
In this section we discuss the phase structures of these theories and clarify when the symmetries 
are broken. It turns out that SU{Nc) QCD with /i/, S0{2Nc) YM and Sp{2Nc) YM should be 
equivalent everywhere in T — ^ parameter space. The equivalence to SU{Nc) QCD with //_b is 



not valid outside the BEC/BCS crossover region of other three theories. In § 3.2.4 we discuss a 
possible cure to this. 

Let us start with S0{2Nc) YM with ^b- A crucial difference from QCD is that there is no 
distinction between 'matter' and 'antimatter' because the gauge group is real. In other words, 
'fundamental' and 'antifundamental' representations are equivalent. For this reason, mesons in 
this theory are not necessarily neutral under U{1)b', one can construct 'baryonic mesons' and 
'antibaryonic mesons' out of two 'quarks' and 'antiquarks', respectively. Because they couple to 
fj,B, as we increase the value of fiB the lightest 'baryonic meson' condenses at some point. Then 
the U{1)b symmetry is broken to Z2 and the equivalence to QCD with /j^b fails. (Note that we 
have used Z4 subgroup of U{1)b for the projection.) 

In order to identify the lightest baryonic meson, let us consider the chiral symmetry breaking in 
this theory. When m = fiB = 0, the Lagrangian ((2j) has the SU{Nf)LX SU{Nf)R x C7(l) b x U{1)a 
symmetry at the classical level at first sight. However, chiral symmetry of the theory is known to 
be enhanced to SU{2Nf). Here U{1)a is explicitly broken by the axial anomaly. One can actually 
rewrite the fermionic part of the Lagrangian ([2]) manifestly invariant under SU{2Nf) using the 
new variable ^ = {i^L,(^2ip*R)^ '■ 

Li = i^^a^D^^, (32) 

where o"^ = {ak,—i) with the Pauli matrices ak {k = 1,2,3). The chiral symmetry SU{2Nf) is 
spontaneously broken down to S0{2Nf) by the formation of the chiral condensate {ipip), leading 
to the 2N^ + Nf - I Nambu-Goldstone bosons living on the coset space SU{2N f) / S0{2N f) 
: neutral pions Ha = 'ip^^Pa'ip, 'baryonic pions' (or simply 'diquark') S5 = ip'^Cj^Qsip and 
'antibaryonic pions' = ip^Cj^Qstp* ■ It is easy to see the fate of these bosons under the 
orbifold projection. The projection to QCD with fiB maps neutral pions to pions in QCD, and 
baryonic and antibaryonic pions are projected away. On the other hand, the projection to QCD 
with fif sends neutral/baryonic/antibaryonic pions to vr^, vr"'" and vr", respectively. Therefore the 
(baryonic) pions in S0{2Nc) YM and those in QCD have the same mass m^. In the same way as 
the vr"'" condensation in QCD with /// at = m,r/2, baryonic pions condense at = m7r/2 (Fig. [5] 
and Fig. [6]) . 

At sufficiently large the one-gluon exchange interaction in the '^■i/'-channel is attractive in the 
color symmetric channel, leading to the condensation of the diquark pairing {ip'^Cj^Qs'tp)- This 
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diquark condensate does not break S0{2Nc) symmetry. This BCS pairing has the same quantum 
numbers and breaks the same U{1)b symmetry as the BEC (S5) at smaU fiB, and there should 
be no phase transition for // > m7r/2 along ^ axis. The phase diagram of this theory is similar to 
that of QCD at finite fj,[, as shown in Fig. [5] and Fig. [6j This is because the condensates in two 
theories are related each other through the orbifold projection, and furthermore, the condensation 
does not break the flavor symmetry used for the projection. 

In the same manner, Sp{2Nc) YM and QCD with fij are equivalent everywhere in T — ^ plane; 
see Fig. [7} (For further details, see [6|.) 

QCD with fiB behaves rather differently, because fiB does not couple to mesons. This does not 
lead to a contradiction, however. Because baryons are much heavier than pions, phenomena char- 
acteristic to QCD with fiB (e.g. formation of hadronic matter) takes place after the equivalence 
is gone due to the U{1)b breakdown in SO{2Nc) and Sp{2Nc) Yang-Mills. 
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Tf 



ii^^C-f^JcQAi') = U(1)s restored 




m7r/2 



Figure 7: Phase diagram of Sp{2Nc) gauge theory at finite fiB- = '4'^C^^JcQa4': where Qa 
{A = 1,2, •• • ,Nf{Nf — l)/2) are antisymmetric Nf x Nf matrices in the flavor space. (Figure 
taken from [6].) 

3.2.4 Breakdown of the equivalence and its cure 

As discussed in [10], the symmetries used in the orbifold projection must not be spontaneously 
broken in order for an orbifold equivalence to hold non-perturbatively. In our case, the Z4 E U{1)b 
of SO{2Nc) and Sp{2Nc) YM, which is necessary to relate them to QCD with fiB, breaks to Z2 
when n > 771^/2 due to 'baryonic pion' condensation. 

In order to recover the equivalence, one can add a deformation term to the S0{2Nc) theory to 
prevent baryonic pion condensation and protect U{1)b- A natural choice for a deformation term 
is 



where A ~ Aqcd- The deformed theory should be viewed as an effective field theory, defined 
with an implicit (lattice) cutoff, since the deformation is an irrelevant operator. At large-A'^c; 
the expectation value of this term factorizes to {Sab) and its complex conjugate, and hence the 
baryonic pion condensation costs C(c^) of energy. Therefore, provided c is large enough, the 
deformation prevents baryonic pion condensation, saving the validity of the orbifold equivalence 
between the deformed S0{2Nc) gauge theory and large Nc QCD. 

Once /i becomes as large as Aqcdi other baryonic mesons might condense, breaking U{1)b- 
Whether this actually happens can only be answered by lattice simulations. If some baryonic 
mesons do condense, one can introduce deformations {e.g., of the form P^^P"-^ where Pab = Qa^Qb 
for the pseudoscalar baryonic mesons) to prevent it, just as was done for baryonic pions above. 
Note that the deformation terms do not survive the orbifold projection to QCD. Thus the value 
of c does not affect the correlation functions of neutral operators so long as it is large enough to 
prevent baryonic pion condensation. 

In order for the numerical simulation, one has to 'integrate in' the auxiliary fields so that the 
action becomes bilinear in the fermion. Whether or not the auxiliary fields introduce the sign 




(33) 



a,b 
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problem is not an easy question. In |5J it has been shown that a similar deformation can be 
introduced without spoiling the sign- free nature. For more on the deformation, see |14j . 



4 Demonstration of the nonperturbative equivalence in solvable 
models 

If the orbifold equivalence holds between the gauge theories, it is natural to expect that the 
equivalence should hold also in solvable toy models which are believed to capture essential features 
of the gauge theories. In this section, we consider the chiral random matrix theory (RMT) |15| . 
and the holographic D3/D7 model which can be studied via AdS/CFT correspondence. We show 
that the perturbative proof applies also in these cases, and furthermore, we confirm the equivalence 
at nonperturbative level. These results strongly suggest that the orbifold equivalence of the gauge 
theories hold nonperturbatively. 



4.1 Chiral random matrix theories 

The partition function of the RMT is given by an integral over a Gaussian random matrix ensem- 
ble, instead of the average over the gauge field of the original Yang-Mills action: 

Nf 

Z = I dci,J]detPe-^^'*=^*'*, (34) 
•' i=i 

where $ is an A/" x (iV + z/) random matrix element, N is the size of the system, v is the topological 
charge, and the Dirac operator T> with quark mass m/ is given by 

. ^ + (35) 

We also introduced a suitable normalization with the parameter G in the Gaussian. Note that 
there is no spacetime coordinate in the theory; the size of the matrix N corresponds to the 
spacetime volume. It is taken to infinity in the end, corresponding to the thermodynamic limit. 
Depending on the anti-unitary symmetries of the Dirac operator, $ is chosen as the real. 



complex, or quaternion real [see (40) for the definition] matrices denoted by the Dyson index 
/3 = 1, /3 = 2, and /3 = 4, respectively. The value of /3 corresponds to the degrees of freedom per 
matrix element. QCD and QCD-like theories corresponding to each universality class are [19]: 

• /3 = 1 : two-color QCD and Sp{2Nc) gauge theory. 

• 13 = 2: SU{N^) QCD {N^ > 3). 

• /3 = 4 : SU{Nc) QCD with adjoint fermions and SO{Nc) gauge theory. 

The effect of temperature T can be incorporated as the (first) Matsubara frequencies by changing 
^ ^ fi + iT for one half of the determinant and ^ ^ jj, — iT for the other half of the determinant 
in the simplest model [18]. 
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The RMT is exactly equivalent to QCD in the e-regime |20] 

-^<.L<.—, fiL<.l, (36) 

where L is the typical scale of the system, because QCD reduces to a theory of zero momentum 
modes of pions. In this regime, the system has a universality, i.e., the dynamics depends only 
on the symmetry breaking pattern and is independent of the microscopic details; QCD can be 
replaced by the RMT with the same global symmetry breaking pattern. Outside the e-regime, the 
universality is lost. However, the RMT is still useful as a schematic model to study the qualitative 
properties of QCD such as the phase structure at finite T and [18]. The advantage of the RMT 
is that it can be solved analytically although QCD itself cannot be. 

As mentioned above, the RMT has the size of the matrix N, which should be taken to infinite 
(thermodynamic limit) in the end. In this sense, the RMT is a "large-A^" matrix model, and 
hence, the perturbative proof of the orbifold equivalence given in Sec. |3] is directly applicable. 
Note that the size of the random matrix is not related to the number of color Nc- The RMT is 
analytically solvable and hence the nonperturbative orbifold equivalence can be checked explicitly. 



4.1.1 Orbifold projections in the chiral random matrix theories 

In this section, we construct the orbifold projections in the chiral random matrix theories (RMTs) 
between /3 = 4, /3 = 2, and (3 = 1. Thereby a class of observables in the RMTs with the different 
Dyson indices are found to be identical to each other. In the following, we will concentrate on the 
RMT at finite n and T = 0, which can be easily generalized to nonzero T. For simplicity, we set 
V = Q and quark masses to be common, mj = m. 

The relationships between these RMTs via orbifold projections are summarized in Fig. |8j We 
start with the /3 = 4 or /3 = 1 RMT at finite with the size of $ being 2N, and define the 
orbifold projection to the /3 = 2 RMT at finite fiB or fii with the size A^. 

The action of the /3 = 4 RMT is given by 

d^d^ e"^, S = Sb + Sf, (37) 
where 

Nf 



and 



V -«>T +/il2Ar ml2N J 

Here are complex Grassmann 4A^-component vectors and $ is a 2A^ x 2A^ quaternion real 
matrix of the form: 



/u=0 
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(3=4 RMT^_____ (3=1 RMT 

with i^B with |iB 




(3=2 RMT^ ^ (3=2 RMT 

withiiB ' ' withiai 



Figure 8: Relations between (3 = 2 RMT at finite fis or /i/ and /3 = 4 and /3 = 1 RMTs at finite 
Hb througli orbifolding. (3 = 2 RMT at small and large fij can be obtained from (3 = 4 and /3 = 1 
RMTs at small and large /u^. (3 = 2 RMT at small fiB can also be obtained from (3 = A RMT at 
small /is, while (3 = 2 RMT at large /xs inside the BEC-BCS crossover region cannot. 

where a'^ are N x N real matrices and o"^ = (— i, a^) with Pauli matrices a^. 
For the bosonic matrix <i>, we impose the projection condition 

J,$J-i = = ( 1^ "^"^ ) • (41) 

Then we obtain 

= ( -I' a" ) • (^^> 

which is equivalent to two copies of a x complex matrix after a unitary transformation 



where 4> = + 10^. The bosonic part of the action is mapped to the one for the (3 = 2 RMT. Note 



that the factor 2 in ( 15 ) is reflected in the difference of normalization between (3 = 4 and (3 = 2 
in (38) if the trace for /3 = 4 is defined as the so-call "QTr" which is one-half the usual trace. 

In order to define a projection for the fermions, we write ^ by using two 2A^-component 
fermions tpn and 1(^1 as 
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Here and ipL are further decomposed into two A^-component fermions 



^ _ f <pl 



where the flavor index is suppressed for simplicity. Then it is straightforward to check that the 
projection 

IpR = UjJclpR, IpL = i^Jc^L, (46) 

where w = e*'^/^ as before, provides us with the /3 = 2 RMT at finite baryon chemical potential. In 
a similar manner, the (3 = 2 RMT at finite isospin chemical potential can be obtained by imposing 
the projection condition 

Jc%^rJ-^ = ^R. Jc^lJ-^ = ^L. (47) 

The /3 = 2 RMT can also be obtained from the /3 = 1 RMT. We start with the action of the 
/3 = 1 RMT given by (38), but this time <1> is a 2N x 2 A'" real matrix, which can be parametrized 
as 

where are N x N real matrices. Note that the only change in this expression compared with 
(40) is that the factors i in front of and are absent. Then one can easily find that the same 
projection conditions for <1> and ^ in the previous subsection gives the (3 = 2 RMT at finite hb 
or finite /i/. 



4.1.2 Explicit demonstration of the nonperturbative equivalence 

The orbifold equivalence in the RMT predicts that the /3 = 4 and (3 = 1 RMTs at finite fiB and 
(3 = 2 RMT at finite fij are equivalent to each other for any m, T and /i in the neutral sector, 
to all order in perturbation theory. Outside the (baryonic) pion condensation phase, the above 
three theories must also be equivalent to the (3 = 2 RMT at finite ^u^. Below we verify the 
nonperturbative equivalence by computing the effective potential of each RMT explicitly. 

For the (3 = 2 and (3 = 1 RMTs, the effective potentials are computed in \W\ I17| . and the 
coincidence of the effective potential of the (3 = 2 RMT at finite fj,B (at /i/ = 0) and that of the 
(3 = 2 RMT at finite fij (at fiB = 0) outside the pion condensed phase has been pointed out. 
Below we first summarize the effective potential of the Nf = 2 RMT with the quark mass m/ 
and the chemical potential fif for each flavor, f = 1,2. The average baryon and isospin chemical 
potentials fis and fij are defined as 

= ^ = 2(/^i+^2), (49) 
= Y = 2(^1-/^2)- (50) 

We denote the parameters of the RMT representing the chiral condensate o"/, pion condensate p, 
and diquark (baryonic pion) condensate A, and their sources as mj, A, and J, respectively. Then 
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we show that the equivalence holds nonperturbatively between RMTs as a natural consequence of 
the orbifold projections. The importance of the unbroken projection symmetry will become clear 
through the argument. 

Effective potential of /3 = 4 

We first consider the /3 = 4 RMT with degenerate quark mass mf = m at finite baryon chemical 
potential /xj = fiB- The effective potential is [6] 



(5=4 



m\ 2 / J 



2 ^ ln[4cj2 + 4A2 - {flB ± iTf]- 



± 



The chiral condensate and the diquark condensate are expressed using a and A as 

1 



4N 
1 

w 



dm In Zj3=4 

dj In Z«=4 



m=0 



J=0 



-4GV;3=4, 



-4G2A^=4. 



(51) 



(52) 
(53) 



Effective potential of /3 = 2 

The effective potential of /3 = 2 is [16] 

G2[(f7i - mi)2 + (as - ms)^ + 2(p - A) 
1 



^13=2 



-\Y^\n[{al+^il± iT) ((72 - T iT) + p^] [(^i - //i ^ iT) (as + /ia ± iT) + 



The chiral condensate and pion condensate are 



mi=0 



(54) 

(55) 
(56) 



A=0 

long as p = (i.e., outside the pion condensed phase), the potential (54) is a function 
h ui)^ and /io = (us — /2/)^. Therefore, 



Note that, as long „„ ^ „ ^ , 

of fif = {fLB + P-i)^ and /ig = (pb — A/)^- Therefore, 



= 0) = n^=2(M/ = PS = 0) for /) = 0. 



(57) 



Here p = means that the projection symmetry, which is used for the orbifolding in Sec. 4.1.1 is 
not broken. From (57), the magnitude of the chiral condensate a and the critical temperature of 
chiral phase transition T'^ in each theory coincide, 



o-/3=2(mb)U/=o = o-/3=2(p/)Ub=o for p = 0, 
Tp=2iP'B)\fij=o = T^=2iP'i)\iiB=o for /' = 0' 



(58) 
(59) 
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as a consequence of the orbifold equivalence. Especially, this shows that the phase-quenched 
approximation for o"/3=2(/^s) and Tp^^iUs) works outside the pion condensed phase, as mentioned 
in US]. 



Note that, even though the effective potentials are identical for p = in (57), the partition 



functions themselves are not generally the same. This is because the pre-exponential factor also 
contributes to the partition function, which is not taken into account in computing the effective 
potential. Therefore, the sign problem measured as the phase of the partition function can be 
severe inside as well as outside the pion condensed phase |21j . The result here shows that the 
phase-quenched approximation is exact for the observables above independently of the severity of 
the sign problem, as long as p = 0. 

Effective potential of /3 = 1 

The effective potential of /3 = 1 is [T7] 



^=1 = G\ai-mif + {a2-m2f + 2{p-Xf + 2{^-jf] 
"i M [(^^1 + Ail ± iT){a2 -fi2T iT) + + A^] 

X [(cJi - ^1 ± iT){a2 + /U2 =F iT) + + A^] + 4AViAi2} 

x{[(ai - ^1 ^ iT){a2 + /i2 ± iT) + p^ + A2] 

X (cji + /ii =F iT){a2 - /X2 ± iT) + p^ + A^] + 4AVi/W2}. 

The chiral condensate, pion condensate, and diquark condensate are 

1 



{uu)p=i 



2N 
1 

w 

1 

In 



dmi In Zi3=i 
dx In Z/3=i 
dj In Zi3=i 



mi=0 



A=0 



= -G cr/3=l, 
/3=1- 



j=o 



The potential (60) has the symmetry 

0/3=i(A, p, pi,p2) = f^/3=i(/0, -A, /ii, -p2), 



(60) 

(61) 
(62) 
(63) 

(64) 



due to the pB l^l symmetry for 13=1. Note that this symmetry has nothing to do with the 
orbifold equivalence. 

Nonperturbative orbifold equivalence between /3 = 4, /3 = 2, and /3 = 1 



By comparing (51), (54), and (60), and by using the pB ^ fJ-i symmetry for /3 = 1, one finds the 



relation (note that A = Oat pB = and /) = at pi = 0): 

17^=4 (2cr/3=4, 2AfS=4)\fig=ij,^fij=o 



2f^/3=2(0"/3=2, P/3=2)U/=/i,/iB=0 
20^=l(cr^=l, A^=l)|^g=^,^j=0- 



(65) 



Unlike the relation (57), this is valid not only for p = (or A = 0) but also for p ^ (or A 7^ 0). 
This is expected because the condensation does not break the projection symmetry, as discussed 
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in § [3j The factor 2 of the effective potentials, mentioned in (15), comes from the fact that the 
/3 = 4 RMT with the size of ^ being 2N twice more degrees of freedom compared to the (3 = 2 
or /3 = 1 RMT with the size N . The origin of the factor 2 for the arguments of the potential 
of the /3 = 4 RMT is the same; fermions in /3 = 4 theory has twice larger degrees of freedom. 



{uu)fs=4 = 2{uu) fj=2 = 2{uu)fs=i. The relation (65) leads to the coincidence of the magnitudes of 



the order parameters (up to the factor 2) and the critical temperatures, 

2cr/3=4(/iB = f^)\/lj=0 = 0-/3=2(/i/ = /^)Us=0 = Cri3=l{UB = /^)Ui=0, (66) 
2A/3=4(/iB = ^)|^^=0 = P/3=2(W = At)lMfl=0 = ^P=l{l^B = ^)Uj=0, (67) 

T^=4ifiB = fJ')\fii=o = T^=2ifj-i = /^)Ub=o = T^=iifj-B = ^J')\^lI=o, (68) 
TpL^ifiB = ^)Uj=o = 7'|=2(W = ^)Ufl=o = T^=i{fiB = At)U/=o, (69) 

which are expected as a consequence of the orbifold equivalence. We note that, the equivalence 
of the neutral order parameters and the critical temperatures should be satisfied in the original 
QCD and QCD-like theories as we claimed in ^ while the effective potentials will not necessarily 
coincide in QCD. The RMT has much less degrees of freedom; the effective potential is a function 
of only the neutral order parameters, and furthermore, all the moments are identical due to the 
orbifold equivalence. As a result, the effective potentials must be identical. In QCD and QCD-like 
theories, the effective potentials depend also on non-neutral observables, and hence the effective 
potentials are not identical in general. 



4.2 Holographic models 

Among other interesting toy models are supersymmetric analogues of large- A'^c QCD which 
have gravity dual description. In [22] the D3/D7 system with chemical potential has been studied 
(Fig. [9]). The starting point is 4d A/" = 4 U(2Nc) super Yang-Mills, which is dual to type IIB 
supergravity in AdS^ x S^. We introduce Nj D7 branes winding on three-cycle of S^. Then open 
strings stretching between D3 and D7 can be regarded as 'quarks' with 'flavor symmetry' U{Nf). 
In the 't Hooft limit Nf/Nc — )• 0, D7-branes behave as probes, and their dynamics is described by 
the Dirac-Born-Infeld (DBI) action on the AdS^ x 5^ background. The isospin chemical potential 
m can be introduced through the boundary condition of the gauge field in the DBI action. This 
theory can be projected to an SO{2Nc) theory with fis through an orientifold projection, and 
further down to U{Nc) theory with fis- However because the orientifold projection does not 
change the local structure of the brane configuration, the equation of motion remains untouched 
as long as the projection symmetry is not broken spontaneously. Therefore, the dynamics of 
mesons determined by the DBI action coincide, and the orbifold equivalence follows. Note that 
this is a 'nonperturbative proof, if we assume the AdS/CFT duality holds nonperturbatively. 

Actually these theories are solved in [23] {U{Nc) with /x^) and [Ml HSj {U{Nc) with /i/), and 
as long as the necessary symmetry is intact the equivalence can be seen explicitly. 



5 Conclusion and outlook 

We have pointed out that QCD and various QCD-like theories with chemical potentials are 
equivalent at large-Ac through the orbifold equivalence, at least to all order in perturbation 
theory. QCD with the isospin chemical potential and SO{2Nc) / Sp{2Nc) Yang-Mills with the 
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U(2Nc), 3 adjoints 
for 2Nf fundamentals 



S0(2Nc), 1 adjoint, 2 symmetries 
fig for 2Nf fundamentals 



U(Nc}, 1 adjoint, 2 symmetries 
Jig for Nf fundamentals 



Figure 9: N = \ Supersymmetric version of the orbifold projection which has a holographic 
realization. 



baryon chemical potential are equivalent everywhere in the T-ji plane, and furthermore, they are 
equivalent to QCD with the baryon chemical potential outside the BEC-BCS crossover region. 
We suggested the equivalence to QCD with the baryon chemical potential may extend inside the 
crossover region by deforming the other theories so that the transition to BEC-BCS phase is 
prevented, while keeping the sign-free nature. 

Our result has immediate implication for the study of the chiral and deconfinement transitions 
in high-T, small-/i region. In this region it is reasonable to assume the Nj/Nc expansion is not 
bad and hence we can expect that the Monte-Carlo results of the QCD with isospin chemical 
potential describe the QCD with the baryon chemical potential with rather good accuracy. (Note 



that the deviation is only 0{{Nf /Nc)"^), as explained in § 3.2.2 ) Furthermore, by using the 
S0{2Nc) theory, one can study three-flavor theory without suffering from the sign problem. In 
a similar manner, from simulation results of the two-color QCD and adjoint QCD, which belong 
to the same universality classes as Sp{2Nc) and SO{2Nc) theories, respectively, one can extract 
information relevant for SU{3) QCD with /j,b- |^ Similar study e.g. phase quenched simulation of 
SU{3) 3-flavor QCD has been performed [26] and the result suggests the QCD critical point does 
not exist. Therefore it is very important to study these sign- free theories numerically, further in 
detail, in order to find (or exclude) the QCD critical point. 
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